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Abstract: A simple evolutionary model for biological ageing is modified 
such that it requires a minimum population for survival, like in human soci- 
ety. This social effect leads to a transition between extinction and survival 
of the species. 

Most computer simulations of biological ageing use the Penna model with 
several free parameters [1,2]. Recently Stauffer [3] in a general review of the 
Penna model suggested a simpler alternative. (For other applications of this 
simpler model see appendix 1). For acceptance in the medical community [4] 
it seems important to include specific human aspects, like the social struc- 
tures of technological societies. Such an attempt is made here, similar in 
spirit but different in details to the simulation of wolf packs in the Penna 
model [5]. 

In the simple alternative [3], each living individual every year produces 
at most one offspring which apart from small mutations by ±1 year has the 
same mimimum reproduction age and the same genetic death age ad as the 
parent. The time unit is called a "year" here, it may correspond four years for 
humans and to days or other time intervals for other species. The probability 
to have one offspring is assumed to be (1 + 0.08)/(a(i(«) — amii) + 0.08), i.e. 
the birth rate is the smaller the longer the reproductive life of the parent 
is: fecundity-survival trade-off [6]. Individuals die if their genetic death age 
Qdi different for different species, is reached; they also die every year with 
a Verhulst probability V^(A^) = N/N^ax due to food and space restrictions. 



where is the actual population at that time, and N^ax ^'^ input parameter, 
often called the carrying capacity (see appendix 2). 

(No individual is allowed to live beyond 32 "years" and no mutations are 
allowed to make > 32 or < a^; similarly we require < a^, < 32.) 

Now we modify this model to include social effects: humans need other 
humans to survive. Thus not only a too high population is bad (as simulated 
by the above Verhulst factor), but also a too low population. We thus take 

V{N) = Nmin/N + N/Nmax {Nmin <N < N^ax) 

where Nmin{^ ^max) is a new parameter which is zero in the old model [3]. 
Fig.l shows this probability for the parameters of our simulations, and Fig. 2 
presents some simulations. The actual population is not in the minimum 
of V{N). If Nmin is not small enough compared with Nmax, the population 
dies out. A fortran program is available from stauffer@thp.uni-koeln.de, a C 
version from janr@icm.edu.pl. Appendix 3 discusses some parameters. 



Verhulst death probability N_min/N + N/N_max used as input for stability limit simulation, N_min = 1450 
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Figure 1: Example of the Verhulst dying probability, making both a too 
small and a too large population unfavourable. 
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Metastability and stability: N_max = 2000 miliion, N_min = 1451 (decay), 1450 (decay), 1449 (stabie) 
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Figure 2: Transition between stable and decaying populations; N„ 
by about a tenth of a percent between these curves. 
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The transition between a surviving and a dying population is quite sharp, 
Fig. 2. If we wait long enough then even a population which first seemed 
stable can become unstable: Fig. 3 shows the rather fiat curve of Fig.2 for 
35 times more years, when quite suddenly it decays. This sudden decay 
comes from the negative feedback of the social term Nmin/N in our Verhulst 
probability: If happens to have a strong downward fluctuation, this dying 
probability Nmin/N becomes larger, thus further decreasing the population. 
The mortality function versus age, and the distributions of and are 
not affected by our age-independent Verhulst factor, and thus the same as in 
Ref. 3. 

(Introduction of a maximum reproduction age can easily kill the popu- 
lation; setting it to 30 when the maximum possible age is 32 reduces the 
average age of genetic death from 18.6 to 18.4.) 

Similar effects are seen for the traditional Penna model [1,2], modified by 
the above Verhulst factor [7] . 



3 



1.2e+06 



N_min = 1449 and N_max = 200 million: Nucleation after long metastable state 



200000 400000 600000 800000 1e+06 1.2e+06 1.4e+06 

time 

Figure 3: Metastable plateau decays after 1.4 million years (continuation of 
the flat curve of Fig. 2). 

In the following appendices we give more details on some of the above pa- 
rameters and assumptions, and summarize other applications of this model. 

We thank T. Klotz for suggesting this work, and S. Moss dc Olivcira for 
communicating her results [7]. JPR was partially supported by the KBN- 
115/E343/2001/ICM grant. 

Appendix 1: Other applications 

Before medical conclusions are drawn from these simulations, more expe- 
rience needs to be gained with this model in other aspects. Thus we shortly 
summarize from the literature other tests of the simple ageing model of ref.[3] . 

The catastrophic senescence of Pacific Salmon, that is the fish dying soon 
after producing offspring, is one of the crucial facts of nature which a theory 
has to explain. It comes out reasonably [8] in the simple ageing model, if one 
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assumes that the fish have offspring only once in their hfes, and all at the 
same age. However, the transition between survival before reproduction, and 
death after reproduction was not as sharp in [8] as in the earlier simulations 
with the more complicated Penna model [9]. 

In 1993 the Northern Cod vanished off the coast of Newfoundland (Grand 
Banks), and the Penna model [10] could explain this effect through over- 
fishing. Fig.4 shows that also the simpler model [3] can explain extinction 
through this overfishing: First we let the fish reach their stationary age dis- 
tribution without any human fishing. Then, after 10000 iterations, humans 
take away one percent of the fish population, resulting in a drop of the sur- 
viving population to a lower but still stable value (center part of Fig.4). 
Finally, after 20000 iterations, another one percent are fished away, and now 
the population shrinks exponentially to zero in the right-most part of this 
figure. 

Vanishing of Northern Cod through overfishing : zero, one and two percent fishing 
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Figure 4: Overfishing in simple model, without social effects. First part: 
no fishing; second part: fishing of one percent gives lower but still stable 
population; third part: fishing increased to two percent leads population to 
extinction. 
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The Penna model was put onto a lattice [11] to simulate e.g. the effect 
that a tree cannot grow on the site where another tree already grows. Similar 
results were obtained [12] if the simple model [3] is put onto a lattice. This 
lattice model allows to avoid an explicit Verhulst factor (see next appendix), 
as pointed out by Makowiec [11]. 

The mortality not only of humans increases roughly exponentially with 
age (Gompertz law) if childhood diseases are overcome. This increase was 
found nicely in the Penna model [1,2] but is followed less accurately in the 
simple model; the latter results follow [12] the Gompertz law better if death 
does not occur abruptly for all individuals at their genetic death age, but 
with some probability increasing in a narrow age interval around the genetic 
death age [13]. Alternatively, reasonable agreement with the Gompertz law 
is obtained by enhancing the birth rate for old individuals [14]. 

Generalization to sexual reproduction also gave reasonable results similar 
to the esexual case [12]. In summary, inspite of its simplicity the model of 
ref.3 passed already the most crucial tests, which have been before apphed 
successfully to the more complicated Penna ageing model. However, some of 
the results are not as good as for the Penna model. 

Appendix 2: Verhulst Factor 

The simple asumption that each individual gives rise to b births in each 
small time interval, and dies with a constant probability d during each time 
interval, leads to a differential equation dN{t)/dt = aN[t) or discrete equa- 
tion N{t + 1) — N{t) = aN{t) for the number N of individuals in a population, 
with a — b — d. This differential equation is solved by an exponential func- 
tion, N{t) = N{0) exp(at), leading to extinction (a < 0) or unlimited growth 
(a > 0). The latter is unrealistic, and thus since the 19th century the equa- 
tion is improved to dN{t)/dt = aN{t)V{t) where V{t) = 1 - N/N^ax is 
the Verhulst factor, taking into account limitations in food and space. The 
carrying capacity N^^ax is in this simple case the infinite-time limit of the 
population: N{t —>■ oo) —>■ Nmax- In more complicated models, like the 
present one of the Penna ageing model, this Verhulst factor V = 1 — N/Nmax 
is only one of several elements, and the infinite-time limit of the population 
no longer has to agree with N^ax', nevertheless is a computationally useful 
and biologically realistic ingredient to take into account the limits of growth. 

Special "chaos" effects may occur [15] in the discrete approach N(t + 1) — 
N(t) = aN{t)V{t) (logistic equation) if the population changes drastically 
from one time step to the other. In our simulations the populations varied 
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smoothly with time and obeyed <^ N^ax thus avoiding these comphca- 
tions. They could appear at higher birth rates, while we assumed one birth 
attempt per time step for each suitable individual. 

In most applications as well as in the present paper, this Verhulst factor 
applies equally to all ages. It may be more realistic biologically to apply it 
only to babies, in which case it acts like a reduction of the birth rate. In 
the Pcnna ageing model, the restriction of the Verhulst factor to babies did 
not give major changes in the equilibrium results [16]; it sometimes gives an 
unrealistic time dependence in how the equilibrium population is approached 
and thus was not applied in the present model. 

Appendix 3: Parameters 

The main advantage of this simple ageing model compared to the more 
reahstic Penna model is that it has few free parameters. We set the birth 
rate to one in the sense that each mature individual in each iteration tries to 
have one offspring, with the probability given in the main text and depending 
on the two characteristic ages. Then the population survives but is close to 
extinction, as is perhaps realistic; thus the population N is far below N^ax- 
With a birth rate (in the above sense) of two, the population was higher 
and no longer close to extinction, but then the increase of mortality with 
increasing age turned out to be less realistic. 

Since the model is close to extinction, the ratio N mini ^ max must be taken 
as very small; in order of magnitude the actual population N is the geometric 
mean of Nmin and Nmax, where the Verhulst death probability Nmin/N + 
N/Nmax has its minimum. We took it such that both survival and extinction 
could be observed by slight changes in Nmin- 

The free parameter 0.08 in the birth probability was introduced originally 
to avoid a wrong divergence if = in the birth probability 1/(0^ — a^)- 
It was then taken such that the population is close to extinction and that 
the age distribution does not extend beyond the maximum age of 32 time 
intervals, traditional for many Penna model studies. This parameter 0.08 
from the simple ageing model without the social effects discussed here was 
then taken over without any change. So the new effects seen here are due to 
our introduction of Nmin and not due to this parameter 0.08. When we omit 
this parameter 0.08 and omit also the restriction of a maximum age of 32, 
we still get a stable population and reasonable results in the simple model, 
also with the present social effects. 
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